An Open-System Quantum Simulator with Trapped Ions 



Julio T. Barreiro*'-^, Markus Miiller*'^'^, Philipp Schindler^, Daniel Nigg^, Thomas Monz^, 

Michael Chwalla^'^, Markus Hennrich^, Christian F. Roos^'^, Peter Zoller^'^ and Rainer Blatt^'^ 
^Institut filr Experimentalphysik, Universitdt Innsbruck, Technikerstr. 25, 6020 Innsbruck, Austria 
^Institut fiir Quantenoptik und Quant eninformation, 
Osterreichische Akademie der Wissenschaften, Technikerstr. 21A, 6020 Innsbruck, Austria 
^Institut fiir Theoretische Physik, Universitdt Innsbruck, Technikerstr. 25, 6020 Innsbruck, Austria 
* These authors contributed equally to this work. 

The control of quantum systems is of fundamental scientific interest and promises 
powerful applications and technologies. Impressive progress has been achieved in iso- 
lating the systems from the environment and coherently controlling their dynamics, 
as demonstrated by the creation and manipulation of entanglement in various physi- 
cal systems. However, for open quantum systems, engineering the dynamics of many 
particles by a controlled coupling to an environment remains largely unexplored. Here 
we report the first realization of a toolbox for simulating an open quantum system 
with up to five qubits. Using a quantum computing architecture with trapped ions, we 
combine multi-qubit gates with optical pumping to implement coherent operations and 
dissipative processes. We illustrate this engineering by the dissipative preparation of 
entangled states, the simulation of coherent many-body spin interactions and the quan- 
tum non-demolition measurement of multi-qubit observables. By adding controlled 
dissipation to coherent operations, this work offers novel prospects for open-system 
quantum simulation and computation. 



Every quantum system is inevitably coupled to its 
surrounding environment. Significant progress has been 
made in isolating systems from their enviroment and co- 
herently controlling the dynamics of several qubits [1-4]. 
These achievements have enabled the realization of high- 
fidelity quantum gates, the implementation of small-scale 
quantum computing and communication devices as well 
as the measurement-based probabilistic preparation of 
entangled states, in atomic [5, 6], photonic [7] and solid- 
state setups [8-10]. In particular, successful demonstra- 
tions of quantum simulators [11, 12], which allow one 
to mimic and study the dynamics of complex quantum 
systems, have been reported [13]. 

In contrast, controlling the more general dynamics of 
open systems amounts to engineering both the Hamilto- 
nian time evolution of the system as well as the coupling 
to the environment. Although open- system dynamics in 
a many-body or multi-qubit system are typically associ- 
ated with decoherence [14-16], the ability to design dissi- 
pation can be a useful resource. For example, controlled 
dissipation allows the preparation of a desired entangled 
state from an arbitrary state [17-19] or an enhanced sen- 
sitivity for precision measurements [20]. In a broader 
context, by combining suitably chosen coherent and dis- 
sipative time steps, one can realize the most general non- 
unitary open-system evolution of a many-particle sys- 
tem. This engineering of the system- environment cou- 
pling generalizes the concept of Hamiltonian quantum 
simulation to open quantum systems. In addition, this 
engineering enables the dissipative preparation and ma- 
nipulation of many-body states and quantum phases [21], 
and also quantum computation based on dissipation [22]. 

Here we provide the first experimental demonstration 



of a complete toolbox, through coherent and dissipative 
manipulations of a multi-qubit system, to control the dy- 
namics of open systems. In a string of trapped ions, 
each ion encoding a qubit, we subdivide the qubits into 
"system" and "environment". The system-environment 
coupling is then engineered through the universal set of 
quantum operations available in ion-trap quantum com- 
puters [23, 24] and a dissipative mechanism based on op- 
tical pumping. 

We first illustrate this engineering by dissipatively 
preparing a Bell state in a 2+1 ion system, such that 
an initially fully mixed state is pumped into a given Bell 
state. Similarly, with 4-hl ions, we also dissipatively pre- 
pare a 4-qubit GHZ-state, which can be regarded as a 
minimal instance of Kitaev's toric code [25]. Besides the 
dissipative elements, we show coherent n-body interac- 
tions by implementing the fundamental building block for 
4-spin interactions. In addition, we demonstrate a read- 
out of n-particle observables in a non-destructive way 
with a quantum-nondemolition (QND) measurement of a 
4-qubit stabilizer operator. Altogether, our work demon- 
strates all essential coherent and dissipative elements for 
controlling general open- system dynamics. 

OPEN-SYSTEM QUANTUM DYNAMICS AND 
BELL-STATE "COOLING" 

The dynamics of an open quantum system S coupled 
to an environment E can be described by the unitary 
transformation psE ^ UpseU\ with psE the joint den- 
sity matrix of the composite system S + E. Thus, the 
reduced density operator of the system will evolve as 
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ps = T^^eU PseU^ ' The time evolution of the system can 
also be described by a completely positive Kraus map 

Ps^£{ps) = ^EkpsEl (1) 

k 

with Ek operation elements satisfying J2k E^Ek = 1 [26]. 
If the system is decoupled from the environment, the gen- 
eral map (1) reduces to ps ^ UspsU\^ with Us the uni- 
tary time evolution operator acting only on the system. 

Control of both coherent and dissipative dynamics 
is then achieved by finding corresponding sequences of 
maps (1) specified by sets of operation elements {Ek} and 
engineering these sequences in the laboratory. In partic- 
ular, for the example of dissipative quantum state prepa- 
ration, pumping to an entangled state reduces to 
implementing appropriate sequences of dissipative maps. 
These maps are chosen to drive the system to the desired 
target state irrespective of its initial state. The result- 
ing dynamics have then the pure state as the unique 
attractor, ps ^ Iii quantum optics and atomic 

physics, the techniques of optical pumping and laser cool- 
ing are successfully used for the dissipative preparation 
of quantum states, although on a single-particle level. 
The engineering of dissipative maps for the preparation 
of entangled states can be seen as a generalization of this 
concept of pumping and cooling in driven dissipative sys- 
tems to a many-particle context. To be concrete, we focus 
on dissipative preparation of stabilizer states, which rep- 
resent a large family of entangled states, including graph 
states and error-correcting codes [27]. 

We start by outlining the concept of Kraus map en- 
gineering for the simplest non-trivial example of "cool- 
ing" a system of two qubits into a Bell state. The 
Hilbert space of two qubits is spanned by the four Bell 
states defined as |<I>±) = ^(|00) ± |11)) and = 

^(|01) ± |10)). Here, |0) and |1) denote the compu- 
tational basis of each qubit, and we use the short-hand 
notation |00) = |0)i|0)2, for example. These maximally 
entangled states are stabilizer states: the Bell state |$^), 
for instance, is said to be stabilized by the two stabi- 
lizer operators Z1Z2 and where X and Z de- 
note the usual Pauli matrices, as it is the only two- 
qubit state being an eigenstate of eigenvalue -hi of these 
two commuting observables, i.e. ^1^2!^+) = |^+) and 
XiX2|^+) = |^+). In fact, each of the four Beh states 
is uniquely determined as an eigenstate with eigenval- 
ues =bl with respect to Z1Z2 and X1X2. The key idea of 
cooling is that we can achieve dissipative dynamics which 
pump the system into a particular Bell state, for example 
Ps ^ |^~)(^~|, by constructing two dissipative maps, 
under which the two qubits are irreversibly transfered 
from the -hi into the -1 eigenspaces of Z1Z2 and X1X2. 

The dissipative maps are engineered with the aid of 
an ancilla "environment" qubit [28, 29] and a quantum 
circuit of coherent and dissipative operations. The form 



and decomposition of these maps into basic operations 
are discussed in Box 1. The cooling dynamics are deter- 
mined by the probability of pumping from the -hi into 
the -1 stabilizer eigenspaces, which can be directly con- 
trolled by varying the parameters in the employed gate 
operations. For pumping with unit probability {p = 1), 
the two qubits reach the target Bell state — regardless of 
their initial state — after only one cooling cycle, i.e., by 
a single application of each of the two maps. In con- 
trast, when the pumping probability is small {p <C 1), 
the process can be regarded as the infinitesimal limit of 
the general map (1). In this case, the system dynam- 
ics under a repeated application of the cooling cycle are 
described by a master equation [30] 

PS = -i\Hs.ps] (2) 
+ ^ {ckPscl - \clckPs - Ps\clc^ . 

Here, Hs is a system Hamiltonian, and are Lind- 
blad operators reflecting the system-environment cou- 
pling. For the purely dissipative maps discussed here, 
Hs = 0. Quantum jumps from the -hi into the -1 
eigenspace of Z1Z2 and X1X2 are mediated by a set of 
two-qubit Lindblad operators (see box 1 for details); here 
the system reaches the target Bell state asymptotically 
after many cooling cycles. 

EXPERIMENTAL BELL-STATE COOLING 

The dissipative preparation of n-particle entangled 
states is realized in a system of n-hl ^^Ca+ ions conflned 
to a string by a linear Paul trap and cooled to the ground 
state of the axial centre-of-mass mode [31]. For each ion, 
the internal electronic Zeeman levels ^5/2(772 = —1/2) 
and Si/2{m = —1/2) encode the logical states |0) and |1) 
of a qubit. For coherent operations, a laser at a wave- 
length of 729 nm excites the quadrupole transition con- 
necting the qubit states (S'1/2 ^ ^5/2)- A broad beam 
of this laser couples to all ions (see Fig. la) and realizes 
the collective single-qubit gate Ux{0) = exp(— i| J^- Xi) 
as well as a M0lmer-S0rensen [32] (MS) entangling oper- 
ation Ux^{0) = ex.p{—i^{Y^.Xi)'^) when using a bichro- 
matic light fleld [33]. Shifting the optical phase of the 
drive field by 7r/2 exchanges Xi by Yi in these operations. 
As a figure of merit of our entangling operation, we can 
prepare 3 (5) qubits in a GHZ state with 98% (95%) fi- 
delity [34]. These collective operations form a universal 
set of gates when used in conjuction with single-qubit ro- 
tations Uzi{0) = exp(— i|Z^), which are realized by an 
off-resonant laser beam that can be adjusted to focus on 
any ion. 

For engineering dissipation, the key element of the 
mapping steps, shown as (i) and (iii) in Box 1, is a sin- 
gle MS operation. The two-qubit gate, step (ii), is real- 



Box 1: Engineering dissipative open-system dynamics 



Dissipative dynamics which cool two qubits from an ar- 
bitrary initial state into the Bell state |^~) are realized 
by two maps that generate pumping from the +1 into the 
-1 eigenspaces of the stabilizer operators Z1Z2 and X1X2: 






























-1 





For Z1Z2, the dissipative map pumping into the -1 
eigenspace is ps ^ ^{ps) — EipsE\ + E2psE\ with 
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The map's action as a uni-directional pumping process can 
be seen as follows. Since the operation element Ei contains 
the projector Z1Z2) onto the +1 eigenspace of Z1Z2, 

the spin flip X2 can then convert +1 into -1 eigenstates of 
Z1Z2, e.g., 1$^) 1-^ 1^^)- In contrast, the -1 eigenspace 
of Z1Z2 is left invariant. In the limit p ^ 1, the repeated 
application of this map reduces the process to a master 
equation with Lindblad operator c — §-^2(1 — Z1Z2). 

We implement the two dissipative maps by quantum 
circuits of three unitary operations (i)-(iii) and a dissipa- 
tive step (iv). Both maps act on the two system qubits S 
and an ancilla which plays the role of the environment E\ 
(i) (ii) (iii) (iv) (i) (ii) (iii) (iv) 



|1> 0*- 
1 •- 









N 








(P| 













(P))l 



^KSh|i) 



""Kx;(p)j/-- 



Cooling Z1Z2 proceeds as follows: 

(i) Information about whether the system is in the +1 
or -1 eigenspace of Z1Z2 is mapped by M{ZiZ2) onto the 
logical states |0) and |1) of the ancilla (initially in |1)). 

(ii) A controlled gate C{p) converts +1 into -1 eigen- 
states by flipping the state of the second qubit with prob- 
ability p, where 

C{p) = |O)(O|o0t/x2(p) + |l)(l|o0l, 

with Ux2 (p) — exp(mX2) and p = sin^ a. 

(iii) The initial mapping is inverted by M~^(ZiZ2). At 
this stage, in general, the ancilla and system qubits are 
entangled. 

(iv) The ancilla is dissipatively reset to which carries 
away entropy to "cool" the two system qubits. 

The second map for cooling into the -1 eigenspace of 
X1X2 is obtained from interchanging the roles of X and 
Z above. 

The engineering of dissipative maps can be readily gen- 
eralized to systems of more qubits. As an example, dissi- 
pative preparation of n-qubit stabilizer states can be real- 
ized by a sequence of n dissipative maps (e.g. for Z1Z2 and 
X1X2X3X4 pumping), which are implemented in analogy 
to the quantum circuits for Bell state cooling discussed 
above: 
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ized by a combination of collective and single-qubit op- 
erations. The dissipative mechanism, step (iv), is here 
carried out on the ancilla qubit by a reinitialization into 
as shown in Fig. lb. Another dissipative process [35] 
can be used to prepare the system qubits in a completely 
mixed state by the transfer |0) (|0) + \S'))/V2 fol- 
lowed by optical pumping of \S') into |1), where \S') is 
the electronic level Si/2{m = 1/2). 

Qubit read-out is accomplished by fluorescence detec- 
tion on the Si/2 ^ P1/2 transition. The ancilla qubit can 
be measured without affecting the system qubits by ap- 
plying hiding pulses that shelve the system qubits in the 
1^5/2 state manifold during fluoresence detection [36]. 

We use these tools to implement up to three Bell-state 
cooling cycles on a string of 2+1 ions. Starting with the 
two system qubits in a completely mixed state, we cool 
towards the Bell state |^~). Each cooling cycle is ac- 
complished with a sequence of 8 entangling operations, 
4 collective unitaries and 6 single-qubit operations; see 
the Supplementary Information. The cooling dynamics 
are probed by quantum state tomography of the system 



qubits after every half cycle. The reconstructed states 
are then used to map the evolution of the Bell-state pop- 
ulations. 

In a first experiment, we set the pumping probability at 
p = 1 to observe deterministic cooling, and we obtain the 
Bell-state populations shown in Fig. 2a. As expected, the 
system reaches the target state after the first cooling cy- 
cle. Regardless of experimental imperfections, the target 
state population is preserved under the repeated appli- 
cation of further cooling cycles and reaches up to 91(1)% 
after 1.5 cycles (ideally 100%). In a second experiment 
towards the simulation of master-equation dynamics, the 
probability is set at p = 0.5 to probe probabilistic cooling 
dynamics. The target state is then approached asymp- 
totically (Fig. 2b). After cooling the system for 3 cycles 
with p = 0.5, up to 73(1)% of the initially mixed popu- 
lation cools into the target state (ideally 88%). In order 
to completely characterize the Bell- state cooling process, 
we also perform a quantum process tomography [26] . As 
an example, the reconstructed process matrix for p = 1 
after 1.5 cycles (Fig. 2c) has a Jamiolkowski process fi- 
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FIG. 1. Experimental tools for the simulation of open 
quantum systems with ions, a, The coherent component 
is reahzed by cohective {Ux ,Uy ^Ux-^ ^Uy^) and single-qubit 
operations (Uzi) on a string of ^"^Ca^ ions which consists of 
the environment qubit (ion 0) and the system qubits (ions 
1 through n). b, The dissipative mechanism on the anciUa 
qubit is reahzed in the two steps shown on the Zeeman-spht 
^°Ca^ levels by (1) a coherent transfer of the population from 
|0) to \S') and (2) an optical pumping to |1) after a transfer 
to the 4^Pi/2 state by a circularly-polarised laser at 397 nm. 
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FIG. 2. Experimental signatures of Bell-state cooling. 

Evolution of the Bell-state populations |$^) (down triangles), 
|$~) (circles), |^^) (squares) and |^~) (up triangles) of an 
initially mixed state under a cooling process with probability 
a, p = 1 or deterministic and b, p = 0.5. Error bars, not 
shown, are smaller than 2% (Icr). c, Reconstructed process 
matrix x (real part), displayed in the Bell-state basis, describ- 
ing the deterministic cooling of the two ions after one and a 
half cycles. The ideal process mapping any input state into 
the state |^~) has as non-zero elements only the four trans- 
parent bars shown. The imaginary elements of x-> ideally all 
zero, have an average magnitude of 0.004 and a maximum of 
0.03. The uncertainties in the elements of process matrix are 
smaller than 0.01 (Icr). 



delity [37] of 87.0(7)% with the ideal dissipative process 
ps ^ which maps an arbitrary state of the 

system into the Bell state |^~). 



FOUR-QUBIT STABILIZER PUMPING 

The engineering of the system-environment coupling, 
as demonstrated by Bell-state cooling above, can be read- 
ily extended to larger n-qubit open quantum systems. 
We illustrate such an engineering experimentally with 
the dissipative preparation of a four-qubit Greenberger- 
Horne-Zeilinger (GHZ) state (|0000) + |1111))/V^. This 
state is uniquely characterized as the simultaneous eigen- 
state of the four stabilizers Z^Z^-, Z^Z-^^ Z-^Z/^ and 
X1X2X3X4, all with eigenvalue +1 (see Fig. 3a). There- 
fore, cooling dynamics into the GHZ state are realized 
by four consecutive dissipative steps, each pumping the 
system into the +1 eigenspaces of the four stabilizers. In 
a system of 4+1 ions, we implement such cooling dynam- 
ics in analogy with the Bell-state cooling sequence. Here, 
however, the circuit decomposition of one cooling cycle 
involves 16 five- ion entangling operations, 20 collective 
unitaries and 34 single-qubit operations; further details 
in the Supplementary Information. 

In order to observe this deterministic cooling process 
into the GHZ state, we begin by preparing the system 
ions in a completely mixed state. The evolution of the 
state of the system after each pumping step is char- 
acterized by quantum state tomography. The recon- 
structed density matrices shown in Fig. 3b for the ini- 
tial and subsequent states arising in each step have a 
fidelity, or state overlap [38] , with the expected states of 
{79(2),89(1),79.7(7),70.0(7),55.8(4)}%; see Supplemen- 
tary Information for further details. Since the final state 
has a fidelity with the target GHZ state greater than 50%, 
the initially mixed state is cooled into a genuinely four- 
particle entangled state [39]. The pumping dynamics is 
clearly reflected by the measured expectation values of 
the stabilizers ZiZ^ {ij = 12,23,34,14) and X1X2X3X4 
at each step, as shown in Fig. 3c. 

Although the simulation of a master equation requires 
small pumping probabilities, as an exploratory study, we 
implement up to five consecutive XiX2X3X4-stabilizer 
pumping steps with two probabilities p = 1 and 0.5, for 
the initial state |1111). The measured expectation val- 
ues of all relevant stabilizers for pumping with p = 1 
are shown in Fig. 3d. After the first step, the stabilizer 
X1X2X3X4 reaches an expectation value of -0.68(1); af- 
ter the second step and up to the fifth step, it is preserved 
at -0.72(1) regardless of experimental imperfections. 

For XiX2X3X4-stabilizer pumping with p = 0.5, the 
four-qubit expectation value increases at each step and 
asymptotically approaches -0.54(1) (ideally -1, fit shown 
in Fig. 3d). A state tomography after each pumping step 
yields fidelities with the expected GHZ-state of {53(1), 
50(1), 49(1), 44(1), 41(1)}%. From the reconstructed 
density matrices we determine that the states generated 
after one to three cycles are genuinely multi-partite en- 
tangled [40]. 
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FIG. 3. Experimental signatures of four-qubit stabilizer pumping, a, Schematic of the four system qubits to be cooled 
into the GHZ state (|0000) + |llll))/v^, which is uniquely characterized as the simultaneous eigenstate with eigenvalue +1 
of the shown stabilizers, b, Reconstructed density matrices (real part) of the initial mixed state pmixed and subsequent states 
pi, 2, 3, 4 after sequentially pumping the stabilizers Z1Z2, ^2^3, ^3^4 and X1X2X3X4. Populations in the initial mixed state 
with qubits i and j antiparallel, or in the -1 eigenspace of the ZiZj stabilizer, disappear after pumping this stabilizer into 
the +1 eigenspace. For example, populations in dark blue dissappear after ZiZ2-stabilizer pumping. A final pumping of the 
stabilizer builds up the coherence between |0000) and |1111), shown as red bars in the density matrix of p4. c, 

Measured expectation values of the relevant stabilizers; ideally, non-zero expectation values have a value of +1. d, Evolution 
of the measured expectation values of the relevant stabilizers for repetitively pumping an initial state |1111) with probability 
p = 0.5 into the -1 eigenspace of the stabilizer The incremental cooling is evident by the red line fitted to the 

pumped stabilizer expectation value. The evolution of the expectation value {X1X2X3X4) for deterministic cooling (p = 1) is 
also shown. The observed decay of (ZiZj) is due to imperfections and detrimental to the pumping process (see Supplementary 
Information). Error bars in c and d, ±lcr. 



COHERENT FOUR-PARTICLE INTERACTIONS 

The coupling of the system to an ancilla particle, as 
used above for the engineering of dissipative dynamics, 
can also be harnessed to mediate effective coherent n- 
body interactions between the system qubits [26, 29]. 
The demonstration of a toolbox for open-system quan- 
tum simulation is thus complemented by adding unitary 
maps ps ^ Us psUl to the dissipative elements described 
above. Here, Us = ex.-p{—irHs) is the unitary time evo- 
lution operator for a time step r, which is generated 
by a system Hamiltonian Hs. In contrast to the re- 
cent achievements [41, 42] of small-scale analog quantum 
simulators based on trapped ions, where two-body spin 
Hamiltonians have been engineered directly [43] , here we 
pursue a gate-based implementation following the con- 
cept of Lloyd's digital quantum simulator [12], where the 
time evolution is decomposed into a sequence of coherent 
(and dissipative) steps. 

In particular, the available gate operations enable an 
experimentally efficient simulation of n-body spin inter- 
actions [44] , which we illustrate by implementing time dy- 
namics of a four-body Hamiltonian Hs = gXiX2XsX4. 



This example is motivated by the efforts to experimen- 
tally realize Kitaev's toric code Hamiltonian [25], which 
is a sum of commuting four-qubit stabilizer operators 
representing four-body spin interactions. This paradig- 
matic model belongs to a whole class of spin systems, 
which have been discussed in the context of topological 
quantum computing [45] and quantum phases exhibiting 
topological order [46]. 



The elementary unitary operation Us can be decom- 
posed into a compact sequence of three coherent opera- 
tions, as explained in Fig. 4a. In an experiment carried 
out with 4+1 ions, we apply Us for different values of r 
to the system ions initially prepared in |1111). We ob- 
served coherent oscillations in the subspace spanned by 
1 0000) and |1111), as shown in Fig. 4b. We character- 
ize our implementation of Us by comparing the expected 
and measured states, determined by quantum state to- 
mography, for each value of r. The fidelity between the 
expected and measured states is on average 85(2)%. 
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FIG. 4. Coherent simulation of 4-body spin interac- 
tions, a, The elementary building block for the simulation 
of coherent evolution Us = ex^p(—irHs) corresponding to the 
four-body Hamiltonian Hs — gXiX2X^X4^ is implemented 
by a circuit of three operations: (i) First, a 5-qubit operation 
M {XiX2X-3,X4) ^ here realized by a single entangling 5- ion MS 
gate [7^2 (7r/2), coherently maps the information, whether the 
four system spins are in the +1(-1) eigenspace of 
onto the internal states |0) and |1) of the ancilla qubit. (ii) 
Due to this mapping, all +1 (-1) eigenstates of XiX^X^X^ ac- 
quire a phase /3/2 (— /3/2) by the single-qubit rotation Uz{/3) 
on the ancilla ion. (iii) After the mapping is inverted, the an- 
cilla qubit returns to its initial state |1) and decouples from 
the four system qubit s, which in turn have evolved according 
to Us' The simulation time step r is related to the phase by 
13 = 2gT. b, Experimentally measured populations in state 
1 0000) (up triangles) and |1111) (circles) as a function of /S for 
a single application of Us to the initial state |1111) of the four 
system qubits (error bars, zblcr). The solid lines show the ideal 
behavior. For comparison, the dashes lines indicate these pop- 
ulations for simultaneous single-qubit (one-body) oscillations, 
each driven by the rotation exp(— 



QND MEASUREMENT OF FOUR-QUBIT 
STABILIZER OPERATORS 

Our toolbox for quantum simulation of open systems is 
extended by the possibility of reading out n-body observ- 
ables in a nondestructive way, which we illustrate here for 
a 4-qubit stabilizer operator X1X2X3X4. As above, we 
first coherently map the information about whether the 
system spins are in the +1(-1) eigenspace of the stabi- 
lizer operator onto the logical states |0) and |1) of the 
ancilla qubit. In contrast to the engineering of coherent 
and dissipative maps above, where this step was followed 
by single-and two-qubit gate operations, here we proceed 
instead by measuring the ancilla qubit. 

Thus, depending on the measurement outcome for the 
ancilla, the system qubits are projected onto the corre- 
sponding eigenspace of the stabilizer: ps ^ Pj^psP+/Nj^ 
{P_psP-/N_) for finding the ancilla in |0) (|1)) with 
the normalization factor A^-j- = Ti{P±psP±)' Here, 
P± = 1(1 ± X1X2X3X4) denote the projectors onto the 
±1 eigenspaces of the stabilizer operator. Note that our 
measurement is QND in the sense that (superposition) 
states within one of the two eigenspaces are not affected 
by the measurement. 

In the experiment with 4+1 ions, we prepare differ- 
ent four-qubit system input states (tomographically char- 



acterize in additional experiments), carry out the QND 
measurement and tomographically determine the result- 
ing system output states. 

To characterize how well the measurement device pre- 
pares a definite state, we use as input |1111), which is a 
non-eigenstate of the stabilizer. In this case, when the 
ancilla qubit is found in |0) or |1) the system qubits are 
prepared in the state (|0000) ± |1111))/V2 by the QND 
measurement. Experimentally we observe this behaviour 
with a quantum state preparation (QSP) fidelity [47] of 
^QSP = 73(1)%. On the other hand, for a stabilizer 
eigenstate, the QND measurement preserves the stabi- 
lizer expectation value. Experimentally, for the input 
state (|0011)-|1100))/V2, we observe a QND fidelity [47] 
of Fqnd = 96.9(6)%. For more details see the Supple- 
mentary Information. 

Our measurement of n-body observables is an essen- 
tial ingredient in quantum error correction and quantum 
computing protocols. In contrast to the open-loop exper- 
iments presented here [28] , this ability also enables an al- 
ternative approach for system-environment engineering: 
The outcome from measurements of the environment can 
be classically processed and used for feedback operations 
on the system. This procedure paves the way to closed- 
loop simulation scenarios in open quantum systems. 

OUTLOOK 

Our experimental demonstration of a toolbox of el- 
ementary building blocks in a system of trapped ions 
should be seen as a first, and conceptual step towards 
the realization of an open quantum system simulator, 
with dynamics governed by the interplay of coherent and 
dissipative evolution. Such a quantum device has appli- 
cations in various fields [13] including condensed-matter 
physics and quantum chemistry, and possibly in mod- 
elling quantum effects in biology [48]. In addition to 
quantum simulation, it enables alternative approaches to 
quantum computing [22]. 

Although the present experiments were performed 
with a linear ion-trap quantum computer architecture, 
the ongoing development of two-dimensional trap ar- 
rays [49] promises scalable implementations of Kitaev's 
toric code [25] and related spin models, as discussed in 
the context of topological quantum computing. Follow- 
ing our original proposal [50] , these ideas can be realized 
with neutral atoms in optical lattices and can be easily 
adapted to other physical platforms ranging from optical, 
atomic and molecular systems to solid-state devices. 
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I. BELL-STATE COOLING 



A. Implemented Kraus maps 



The Bell state |^~) is not only uniquely determined 
as the simultaneous eigenstate with eigenvalue -1 of the 
two stabilizer operators X1X2 and Z1Z2 (as mentioned 
in the text), but also by X1X2 and Y1Y2. In the exper- 
iment, we implemented cooling into |^~) by engineer- 
ing the two Kraus maps ps ^ EipsE\ + E2PSE2 and 



B. Circuit decomposition 

The map for pumping into the 
was implemented by the unitary 



-1 eigenspace of X1X2 



X^{7r/2) 



(1) 



(corresponding to steps (i) - (iii) in Box 1) followed by 
an optical pumping of the ancilla qubit to Here, the 
two-qubit controlled gate is 

C{p) = |0){0|o (8) exp(mZi) -h |l){l|o ^ 1 



exp 



-(l + Zo)iceZi 



= UzA-o^) UY{n/2) U^^i'\-a) C/y(-7r/2) (2) 

where U^i^\—a) = exp(i(a/2)XoXi) denotes an MS 
gate acting only on the ancilla and the first system 
qubit. This two-qubit MS gate operation was imple- 
mented in the experiment by the use of refocusing tech- 
niques (Nebendahl et a/., 2009). In more detail, the 
gate ^-^j^2 

was realized by interspersing two of the avail- 
able three-qubit MS gate operations with single-ion light 



2 



shifts on the second system qubit which induces a 7r-phase 
shift between the qubit states. Alternatively, this refo- 
cusing could be avoided, and the sequences further sim- 
plified, by hiding the population of individual ions (here 
the second system ion) which are not supposed to partic- 
ipate in collective coherent operations in electronic levels 
decoupled from the driving laser excitation. More details 
on how to systematically decompose Kraus maps into the 
experimentally available ion-trap gate operations, in par- 
ticular the multi-ion MS entangling gate, can be found 
in (Miiller et al, 2010). 

The circuit decompositions for the experimental imple- 
mentation of the two maps are shown in Fig. SI. We note 
that the circuits have been simplified at the expense of 
implementing in addition in each dissipative map a fiip 
operation Y1I2 on the two system qubits. However, as 
this additional unitary corresponds to one of the stabi- 
lizers into whose -1 eigenspace cooling is performed, this 
does not interfere with the cooling dynamics. 

Cooling with unit pumping probability p = 1 corre- 
sponds to a = 7r/2, whereas p = 0.5 is realized with by 
setting a = 7r/4. In the experiment, the "fundamental" 
MS gate was calibrated to implement [/x2(a/2). The 
fully entangling operation Ux'^ {^l"^) ^t the beginning and 
the end of the sequence Fig. Sla was then implemented 
by applying the Ux'^{ol/2) operation twice (for p = 1) 
or four times (for p = 0.5). The fully entangling oper- 
ations [/y2(7r/2) in Fig. Sib were implemented by two- 
and four- fold application of the "fundamental" MS gate 
with a shifted optical phase of the driving laser (cf. Sec- 
tion 2 in the main text). 
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FIG. SI Experimental sequences for Bell-state cool- 
ing. Pumping into the eigenspaces of eigenvalue -1 of X1X2 
(circuit a,) and Y1Y2 (circuit b,) occurs with a probability p 
in each step, where sin^ a = p. 



C. Additional data 

The initially mixed state was prepared with a fidelity 
of F=99.6(3)% with respect to the ideal state |l4x4- 

Physical process matrices were reconstructed with 
maximum likelihood techniques (Jezek et a/., 2003). An 



error analysis was carried out via Monte Carlo simu- 
lations over the multinomially distributed measurement 
outcomes of the state and process tomography. For each 
process and state, 200 Monte Carlo samples were gener- 
ated and reconstructed via maximum-likelihood estima- 
tion. 



II. FOUR-QUBIT STABILIZER PUMPING 

Expectation values of the stabilizer operators ^1^2 5 
Z2Z3, Z3Z4 and X1X2X3X4 were not determined from 
the reconstructed density matrices of the system qubits. 
Instead, we performed fluorescence measurements in the 
X and Z basis on 5250 copies of the corresponding quan- 
tum states (for p = 0.5 cooling, 2100 copies were mea- 
sured). The error bars were then determined from the 
multinomially distributed raw data. 



A. Cooling 

Cooling into the GHZ state (|0000) + |1111))/V2 was 
realized by a pumping cycle where the four system qubits 
were deterministically pumped into the +1 eigenspaces of 
the stabilizers Z1Z2, ^2^3, Z^Z^ and X1X2X3X4. 

Pumping into the -1 eigenspace of Z1Z2 in the first 
cooling step could be achieved in complete analogy with 
Bell state cooling, i.e. by implementing a dissipative map, 
which only involves operations on the ancilla qubit and 
the system qubits #1 and #2, whereas the system qubits 
#3 and #4 remain completely unaffected. This could 
either be achieved through refocusing techniques or by 
hiding system ions #3 and #4 in electronically decoupled 
states for the duration of the dissipative circuit. 

In the experiment, however, we used a few simplica- 
tions, which are schematically shown in Fig. S2 and listed 
below: 

• For deterministic cooling (p = 1), the inverse map- 
ping step (shown in Box 1) is not necessary and has 
been taken out. 

• In the coherent mapping step (shown in Box 1) the 
information about whether the system ions are in 
a ±1 eigenstate of Z1Z2 is mapped onto the logical 
states of the ancilla qubit. This step ideally only 
involves the ancilla and the system qubits #1 and 
#2. One way to achieve this three-qubit operation 
without affecting the system qubits #3 and #4, is 
to combine the available five-ion MS gate with ap- 
propriately chosen refocusing pulses, i.e. light shift 
operations on individual ions. Those would have to 
be chosen such that ions #0, #1 and #2 become 
decoupled from ions #3 and #4, and furthermore 
residual interactions between ions #3 and #4 can- 
cel out. However, it turns out that residual inter- 
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FIG. S2 Pumping into the +1 eigenspace of the Z1Z2 
stabilizer operator a, Ideally, only the ancilla qubit and the 
two system qubits #1 and #2 are involved in the circuit, b, 
An entangling gate acting on these three ions can be achieved 
by a refocusing technique, where ions #3 and #4 decouple 
from the dynamics. However, the latter ions still become en- 
tangled. However, these residual interactions are not harmful 
to the cooling, as they do not affect the expectation values 
of the other two-qubit stabilizer operators, c, Dashed opera- 
tions in the quantum circuit indicate such residual entangling 
operations. 



actions between ions #3 and #4 can be tolerated: 
although not required for the ZiZ2-pumping dy- 
namics, they are not harmful, as they do not al- 
ter the expectation values of the other two-qubit 
stabilizers Z2Z3 and Z3Z4. In our experiment the 
decoupling of ions #0, #1 and #2 from the ions 
#3 and #4 was achieved by the circuit shown in 
Fig. S2b. 

The additional interactions in the pumping of the 
two-qubit stabilizer operators ZiZj affect the state 
of the system qubits with respect to the four-qubit 
stabilizer X1X2X3X4. However, this effect is not 
detrimental to the cooling, provided the pumping 
into the eigenspace of X1X2X3X4 is performed as 
the final step in the cooling cycle. 

• In the employed sequence, the number of single- 
qubit rotations was reduced wherever possible. Es- 
sential single-qubit light shift operations, such as 
those needed for re- focusing operations, were kept. 

• Local rotations of the system ions at the end of 
a cooling step, which would be compensated at 



the beginning of the subsequent cooling step, were 
omitted when several dissipative maps were applied 
in a row. The corresponding gates of the sequences 
are displayed in blue in Steps 1-3. 

These simplifications allowed us to significantly reduce 
the length and complexity of the employed gate sequences 
for one stabilizer pumping step. The compressed gate 
sequences as used in the experiment are explicitly given 
below. 

Step 1 (pumping into the +1 eigenspace of Z1Z2): 

Uy{-7t/2)UzA-^/2) 
UxiTr/2)Uz,i-n/2)Uxi-7T/2) 
Uz, (tt) Ux2 (7r/4) Uz, (tt) Uz, (tt) Ux^ (7r/4) 
Ux{-T^l2) UzA-^/2) f/zo(-V2) Ux{^l2) 
Ux2 (7r/4) Uz, (tt) Uz, (tt) Ux- (7r/4) 
Uy{7t/2) Ux{-7t/2) Uzo{-7r/2) Ux{7t/2) 

Step 2 (pumping into the +1 eigenspace of Z2Z3): 

Uy{-7t/2)UzA-^/2) 
Ux{TT/2)Uz,{-n/2)Ux{-n/2) 
Uz, (tt) Ux2 (7r/4) Uz, (tt) Uz, (tt) Ux^ (7r/4) 
Ux{-7r/2) Uz,{-^/2) f/zo(-V2) Ux{^l2) 
Ux- (7r/4) Uz, (tt) Uz, (tt) Ux- (7r/4) 
Uy{^/2) Uxi-7T/2) Uzoi-Tr/2) Ux{7t/2) 

Step 3 (pumping into the +1 eigenspace of Z^Z^): 

C/y(-7r/2)f/2,(-7r/2) 
Ux{tt/2)Uz,{-tt/2)Ux{-tt/2) 
Uz, (tt) Ux- (7r/4) Uz, (tt) Uz, (tt) Ux- (7r/4) 
Ux{-T^l2) VzA-T^m f/zo(-V2) C/x(7r/2) 
Ux2 (7r/4) Uz, (tt) Uz, (tt) Ux^ (7r/4) 
C/y (7r/2) f/x(-7r/2) Uz,{-t^I2) Ux{^l2) 



Step 4 (pumping into the +1 
X1X2X3X4): 



eigenspace 



of 



UzA-^/2)Ux{^/2)UzA-^/2) 
Ux- (^/4) Uz, (^) Uz, (^) Ux^ (^/4) 
UzA-^m Ux{-7t/2) Uz,{-7r/2) Ux{n/2) 
Ux2{7T/4)Ux2{7r/4) 

Figure SIO shows the reconstructed density matrices 
(real and imaginary parts) for every step of the cooling 
cycle. The complete circuit decomposition of one cooling 
cycle involves 16 five-ion entangling operations, 28 (20) 
collective unitaries and 36 (34) single-qubit operations 
with (without) optional operations in blue. The reset 
operation involves further pulses not accounted for above. 
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B. Repeated four-qubit stabilizer pumping 

To study the robustness of the dissipative operation, 
we prepared the initial state |1111) and subsequently ap- 
plied repeatedly the dissipative map for pumping into the 
+1 eigenspace of the four-qubit stabilizer X1X2X3X4. 
We observed that after a single dissipative step a non- 
zero expectation value of X1X2X3X4 built up and stayed 
constant under subsequent applications of this dissipative 
map. However, due to imperfections in the gate opera- 
tions, the expectation values of the two-qubit stabiliz- 
ers decreased, ideally they should not be affected by the 
XiX2X3X4-pumping step (see Fig. S3). Interestingly, 
the expectation values of Z\Z^ and decayed signif- 

icantly faster than those for ZxZ^ and Z^Z-^. This decay 
can be explained by the fact that in the gate sequence 
used for pumping into the +1 eigenspace of X1X2X3X4, 
step 4 above, single- ion light-shift operations are applied 
only to the fourth system qubit and the ancilla. This 
indicates that errors in the single-qubit gates applied to 
the fourth system ion accumulate under the repeated ap- 
plication of the dissipative step, and thus affect the sta- 
bilizers Z1Z4 and Z3Z4 which involve this system qubit 
more strongly than the others. This destructive effect 
can be minimized by alternating the roles of the system 
qubit s. 




|1111> Wr Wr Wr mr mr 



FIG. S3 Measured expectation value of stabilizers 
for repeated pumping without sequence optimization. 

The expectation values of Z1Z4 and Z3Z4 show a significantly 
faster decay than those for Z1Z2 and Z2Z3. In every step of 
the cooling, most single-ion light-shift operations are applied 
to the fourth system qubit. 

Such optimization has been done for the dissipative 
dynamics shown in Fig. S4. Here, starting from the initial 
state 1 1111), repeated pumping into the -1 eigenspace of 
X1X2X3X4 has been implemented by the sequence 

Ux^ (^/8) Ux2 {tt/S) Ux2 {tt/S) Ux- (^/8) 

Ux{-^m 

UzA-^/2xp)Ux{^/2)UzM 
[/y2(7r/4 X p) UzoM UzA^) t/y2(7r/4 x p) 
Uy{^/2)Uz,{-^/2)Uy{-^/2) 
Ux- (^/8) Ux^ (^/8) Ux^ (^/8) Ux^ (^/8) 
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FIG. S4 Measured expectation value of stabilizers for 
repeated pumping with sequence optimization. All 

two-qubit stabilizers decay at the same rate during cooling. 
In step 1,2,3,4, and 5 the single-qubit light-shift operations 
were applied on the system qubits 4,3,2,1, and 1, respectively. 

Here, we observed that indeed the expectation values 
of all two-qubit stabilizers decreased at the same pace 
and at a slightly slower rate (see Fig. S4). Upon repeat- 
ing the sequence above 1,2,3,4, and 5 times, we changed 
the operations shown in red to act on qubits 4,3,2,1, and 
1, respectively. The stabilizer expectation values for de- 
terministic cooling, or p = 1, are shown in Fig. S4. 



C. Pushing "anyons" around 

In Kitaev's toric code (Kitaev, 2003), spins are located 
on the edges of a two-dimensional square lattice. The 
Hamiltonian 

H = -g{J2Ap + J2Bp) (3) 

P V 

is a sum of mutually commuting four-qubit stabilizers 
= Iliep^i = Iliev^i^ which describe four- 

spin interactions between spins located around plaque- 
ttes p and vertices v of the lattice. The ground state of 
the Hamiltonian is the simultaneous +1 eigenstate of all 
stabilizer operators. The model supports two types of 
excitations that obey anyonic statistics under exchange 
(braiding), and they correspond to -1 eigenstates of either 
plaquette or vertex stabilizers. 

For a minimal instance of this model, represented by 
a single plaquette of four spins located on the edges, the 
Hamiltonian contains a single four-qubit interaction term 
X1X2X3X4 and pairwise two-spin interactions ZiZj of 
spins sharing a corner of the plaquette. The ground state 
as the simultaneous +1 eigenstate of these stabilizers is 
the GHZ-state (|0000) + |1111))/V^. States correspond- 
ing to -1 eigenvalues of a two-qubit stabilizer ZiZj can be 
interpreted as a configuration with an excitation located 
at the corner between the two spins i and j. Similarly, a 
four-qubit state with an eigenvalue of -1 with respect to 
X1X2X3X4, would correspond to an anyonic excitation 
located at the center of the plaquette. 
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In the experiment we prepared an initial state |0111) 
and then performed the coohng cycle of four determinis- 
tic pumping steps into the +1 eigenspaces of Z1Z2, ^2^3, 
Z3Z4 and X1X2X3X4, using the sequences for Steps 1 to 
4 given in section 11. A. The expectation values of the sta- 
bilizer operators for the initial state and the four spins 
after each pumping step are shown in Fig. S5. The dis- 
sipative dynamics can be visualized as follows: For the 
initial state with (Z1Z2) = — 1 and (Z1Z4) = — 1 a pair 
of excitations is located on the upper left and right cor- 
ners of the plaquette, whereas (X1X2X3X4) = implies 
an anyon of the other type is present at the center of 
the plaquette with a probability 50%. In the first cool- 
ing step, where the first two spins are pumped into the 
+1 eigenspace of Z1Z2, the anyon at the upper right cor- 
ner is dissipatively pushed to the lower right corner of 
the plaquette. In the third step of pumping into the +1 
eigenspace of Z3Z4, the two excitations located on the 
upper and lower lefts corners fuse and disappear from 
the system. In the final step of pumping into the +1 
eigenspace of X1X2X3X4, the anyon with a probability 
of 50% at the center of the plaquette is pushed out from 
the plaquette. 

However, we'd like to stress that borrowing concepts 
from topological spin models, such as anyonic excitations, 
here is merely a convenient language to phrase and visual- 
ize the dissipative dynamics. In the present work with up 
to five ions, we do not explore the physics of topological 
spin models, since (i) in a minimal system of four spins 
the concepts developed for larger lattice models become 
questionable, and more importantly, (ii) during the im- 
plemented cooling dynamics the underlying (four-body) 
Hamiltonian of the model was not present. We rather 
demonstrate the basic tools which will allow one to ex- 
plore this physics once larger, two-dimensional systems 
become available in the laboratory. 

We note that photon experiments have reported the 
observation of correlations compatible with the manipu- 
lations of "anyons" in a setup representing two plaque- 
ttes (Lu et al.^ 2009; Pachos et al.^ 2009). Such exper- 
iments are based on postselection of measurements (as 
in teleportation by Bouwmeester et al.^ 1997), which 
should be contrasted to our deterministic implementation 
of open system dynamics to prepare and manipulate the 
corresponding quantum state (as in deterministic tele- 
portation by Barrett et al, 2004; Riebe et al, 2004). 



I I |U I rII 

^ I I .1 M (X1W4) 

10111) VWlr Pi P2 VW\r P3 Mr P4 

Z/, Z,Z3 Z3Z, X/,X3X, 

FIG. S5 Pushing "anyons" around by dissipation. Mea- 
sured expectation values of stabilizer operators for cooling dy- 
namics of pumping into the +1 eigenspaces of Z1Z2, Z2Z3, 
Z3Z4 and X1X2X3X4, starting in the state |0111). 

ing into the -1 eigenspace of X1X2X3X4. In the context 
of Kitaev's toric code, this state would correspond to 
an excited state. However, as above, we point out that 
the underlying Hamiltonian was not implemented in the 
cooling dynamics. 

The measured expectation values of the stabilizers are 
shown in Fig. S6. The final density matrix, as deter- 
mined from quantum state tomography after the four 
cooling steps, is shown in Fig. S7. This cooling cycle was 
implemented with the same sequences as given for Step 
1 to 4 in section H.A, with the only difference that the 
sign of the phase shift operations displayed in red was 
changed in Steps 2, 3, and 4. This allowed us to invert 
the pumping direction from the +1 into -1 eigenspaces of 
Z2Z2', ^3^4 ^iid X1X2X3X4. 
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FIG. 86 Cooling into an "excited" state. Measured ex- 
pectation values of two- and four-qubit stabilizer operators for 
pumping into the state (|0010) — |1101))/v^, starting from an 
initially four-qubit mixed state. 



D. Pumping into "excited" states 

Starting from an initially fully mixed state of four 
qubits, we also implemented cooling into a different GHZ- 
type state, (|0010) — 1 1101)) /a/2, by a sequence of four 
dissipative steps: 1) pumping into the +1 eigenspace of 
Z1Z2, 2) pumping into the -1 eigenspace of 3) 
pumping into the -1 eigenspace of Z3Z4 and 4) pump- 



Ill. QND MEASUREMENT OF A FOUR-QUBIT 
STABILIZER 

A. Further details 

As shown in Fig. S8, the QND measurement involves a 
mapping step where the information about whether the 
system described by an input density matrix p^^ is in 



6 




FIG. S7 Reconstructed density matrix after the full 
cooling cycle for dissipative preparation of the state 

(|0010) - |1101))/V2. This final state has a fidelity of 
60(2)% with the expected state. This fidelity was determined 
from parity and coherence measurements and analysed with 
bayesian inference techniques as done in (Monz et aL, 2010). 



the system qubits being in a state belonging to the +1 
eigenspace of the stabilizer operator, the ancilla is flipped 
from |1) to |0), whereas it remains in its initial state |1) 
otherwise. 

Subsequently, the ancilla as well as the four system 
qubits were measured. This was done by measuring the 
five ions simultaneously. Alternatively, we first hid the 
four system qubits in electronic levels decoupled from the 
laser excitation, performed the fluorescence measurement 
of the ancilla qubit, then recovered the state of the sys- 
tem qubits and tomographically measured the state of 
the four system qubits. The second approach, where the 
state of the system is not affected by the measurement of 
the ancilla, is of importance if the information from the 
ancilla measurement is to be used for feedback operations 
on the state of the system. 



the +1 / -1 eigenspace of A = X1X2X3X4 is coherently 
mapped onto the internal states |0) and |1) of the ancilla 
qubit, which is initially prepared in |1). Subsequently 
the ancilla qubit is measured in its computational basis, 
leaving the system qubits in a corresponding output state 
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FIG. 88 QND measurement of the four-qubit stabi- 
lizer operator X1X2X3X4. After the coherent mapping 
M{XiX2X3X4), the ancilla qubit is measured. This mea- 
surement was performed both with and without applying ad- 
ditional pulses to hide the populations of the system qubits 
in electronically uncoupled states for the duration of the flu- 
orescence measurement on the ancilla. 

The coherent mapping M{XiX2XsX4) was realized by 
the sequence 

Ux{7T/4)Uzo{7T)Ux{-7r/4) 

Ux^ (^/4) Ux^ (^/4) Uz, (-^/2) Ux^ (^/4) Ux^ (^/4) 
Uy{-^/A)Uz,{^)Uy{^/A) 

which implements 

M(XiX2X3X4) = - ^(Xo + Fo) P+ 
V2 



'V2 



(l-iZo)0P-, (4) 



with P± = ^(111X1X2X3X4) the projectors onto the ±1 
eigenspaces of X1X2X3X4. Equation (4) shows that for 



B. Quantitative analysis of the performance 

To characterize the performance of a QND measure- 
ment for a (multi-) qubit system, a set of requirements 
and corresponding fidelity measures have been discussed 
in the literature (Ralph et al.^ 2006). 

(1) First of all, the measurement outcomes for the an- 
cilla qubit should agree with those that one would expect 
from a direct measurement of the observable A on the in- 
put density matrix. This property can be quantified by 
the measurement fidelity^ 



^M = (v^^+V^^)' 



(5) 



which measures the correlations of the distribution of 
measurement outcomes = of the ancilla 

qubit with the expected distribution p^^ = {p^^^p^I^} di- 
rectly obtained from p^", where p^^ = Tr{|(l ± A)p^^}. 

(2) The QND character, reflected by the fact that the 
observable A to be measured should not be disturbed 
by the measurement itself, becomes manifest in ideally 
identical probability distributions p^^ and which are 
determined from the input and output density matrices. 
These correlations are quantified by the QND fidelity 



^QND 



^in^out 



pin^out 



(6) 



where p^^* = Tr{|(l ± A)p°^*}. 

(3) Finally, by measuring the ancilla qubit the sys- 
tem qubits should be projected onto the corresponding 
eigenspace of the measured observable A. Thus the qual- 
ity of the QND measurement as a quantum state prepa- 
ration (QSP) device is determined by the correlations 
between the ancilla measurement outcomes and the cor- 
responding system output density matrices. It can be 
described by the QSP fidelity 



QSP 



. m out 
"^+^|0),+ 



m out 



(7) 



where p^^Ji^ ^ denotes the conditional probabihty of find- 
ing the system qubits in the +1 (-1) eigenspace of A, pro- 
vided the anciha qubit has been previously measured in 

|o) (|i)). 

The probability distributions for the system in- 
put and output states, the ancilla measurement out- 
come distributions, and the resulting fidelity values 
are summarized in Tables I to IV. The input states 
had a fidelity (Jozsa, 1994) with the ideal states 
(|0000) + |1111))/V2, (|0000) - |1111))/a/2 and (|0011) - 
|1100))/\/2 of 75.3(9), 77.3(8), 93.2(4)%. 

We observe that we obtain higher values for the mea- 
surement and QND fidelities than for the QSP fidelities. 
The latter is relevant in the context of quantum error cor- 
rection or closed- loop simulation protocols or more gener- 
ally whenever the information from the ancilla measure- 
ment is used for further processing of the system output 
state. 

With the additional hiding and unhiding pulses before 
and after the measurement of the ancilla we observe a 
loss of fidelity of a few percent in the QSP fidelities. 
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TABLE I QND probability distributions. Obtained from measurements with hiding of the system ions during the measurement 
of the ancilla. 



input state 


eigenspace 


m m in in in out out out 
Pin Pout P Pni=0 Pni = l P Pni=0 Pni=l 


|0000> + |1111> 


+1 

-1 


0.959(1) 0.847(3) 0.817(9) 0.822(9) 0.618(34) 0.689(12) 0.736(12) 0.359(34) 
0.041(1) 0.153(3) 0.183(9) 0.178(9) 0.382(34) 0.311(12) 0.264(12) 0.641(34) 


lOOOO) - |1111) 


+1 
-1 


0.955(1) 0.169(3) 0.191(10) 0.187(9) 0.328(36) 0.310(11) 0.640(26) 0.242(12) 
0.045(1) 0.831(3) 0.809(10) 0.813(9) 0.672(36) 0.690(11) 0.360(26) 0.758(12) 


|0011) - |1100) 


+1 
-1 


0.978(1) 0.103(2) 0.041(4) 0.035(4) 0.412(47) 0.137(9) 0.476(36) 0.097(7) 
0.022(1) 0.897(2) 0.959(4) 0.965(4) 0.588(47) 0.863(9) 0.524(36) 0.903(7) 



TABLE II QND probabihty distributions. Obtained from measurements without hiding of the system ions during the 
measurement of the anciha. 

input state eigenspace pZt p""""^ pZlo P^Li 



|0000) + |1111) 


+1 
-1 


0.850(3) 0.713(11) 0.789(11) 0.336(30) 
0.150(3) 0.287(11) 0.211(11) 0.664(30) 


|0000> - 


11111) 


+1 

-1 


0.188(3) 0.265(12) 0.504(28) 0.220(11) 
0.812(3) 0.735(12) 0.496(28) 0.780(11) 


looii) - 


|1100) 


+1 
-1 


0.099(2) 0.073(7) 0.416(35) 0.038(5) 
0.901(2) 0.927(7) 0.584(35) 0.962(5) 



TABLE III QND figures of merit. Determined from measurements with hiding of the system ions during the measurement of 
the anciha. Since the state |0011) — |1100) is particularly robust against decoherence, the fidelity Fqsp is higher, as shown for 
8 ions in (Monz et al, 2010). 



input state 


eigenspace 


pin pOut 






1^QND = - 


Fmip'^'^p'^ 


) ^QNr>{P ,p 


) FQSP(l^"^,P?rND) 


lOOOO) + |1111) 


+1 


0.82(1) 0.69(1) 


0.85 


0.74(1) 




0.998(1) 


0.978(5) 


0.72(1) 




-1 


0.18(1) 0.31(1) 


0.15 




0.64(3) 








|0000) - |1111) 


+1 


0.19(1) 0.31(1) 


0.17 


0.64(3) 




0.999(1) 


0.980(5) 


0.74(1) 




-1 


0.81(1) 0.69(1) 


0.83 




0.76(1) 








|0011) - |1100) 


+1 


0.041(4) 0.14(1) 


0.10 


0.48(4) 




0.985(3) 


0.969(6) 


0.86(1) 




-1 


0.959(4) 0.86(1) 


0.90 




0.90(1) 








11111) 


+1 


0.5 0.47(1) 0.50049 


0.70(1) 




1 


0.9992(6) 


0.73(1) 




-1 


0.5 0.53(1) 0.49951 




0.76(1) 
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TABLE IV QND figures of merit. Determined from measurements without hiding of the system ions during the measurement 
of the ancilla. Since the state |0011} — |1100} is particularly robust against decoherence, the fidelity Fqsp is higher, as shown 
for 8 ions in (Monz et al, 2010). 



input state eij 


^enspace 


pOut 




Pqnd=+ 


Pqnd=- 


i^M(p^^p"^ 


) FQND(p^^p°"^ 


) FQSP(l^"^,P?rND) 


lOOOO) + lllll) 


+1 

-1 


0.71(1) 
0.29(1) 


0.85 
0.15 


0.79(1) 


0.66(3) 


0.998(1) 


0.984(4) 


0.77(1) 


lOOOO) - lllll) 


+ 1 
-1 


0.26(1) 
0.74(1) 


0.19 
0.81 


0.50(3) 


0.78(1) 


1.0000(1) 


0.992(3) 


0.73(1) 


lOOll) - |1100) 


+1 
-1 


0.07(1) 
0.93(1) 


0.10 
0.90 


0.42(3) 


0.96(1) 


0.986(2) 


0.996(2) 


0.91(1) 


lllll) 


+1 
-1 


0.52(1) 
0.48(1) 


0.5078 
0.4922 


0.75(1) 


0.73(1) 


0.99994 


0.9996(5) 


0.74(1) 



10 



Bell-State pumping with p = 1 after 1 pumping cycle 




Bell-state pumping with p = 0.5 after 3 pumping cycles 




FIG. S9 Reconstructed process matrices of experimental Bell-state cooling. The reconstructed process matrix for 
p = 1 after 1 (1.5) cycles has a Jamiolkowski process fidelity (Gilchrist et a/., 2005) of 83.4(7)% (87.0(7)%) with the ideal 
dissipative process ps i-^ |^~)(^~| which maps an arbitrary state of the system into the Bell state |^~). This ideal process 
has as non-zero elements only the four transparent bars shown. The reconstructed process matrix for p = 0.5 after 3 cycles has 
a Jamiolkowski process fidelity of 60(1)% with the ideal process Xideai shown [Im(xideai) = 0]. 
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FIG. SIO Ideal and reconstructed density matrices of plaquette cooling. An inital mixed state pmixed is sequentially 
pumped by the stabilizers Z1Z2, ^2^3, ^3^4 and driving the system into the states pi, 2,3,4- 




FIG. Sll Pushing "anyons". Cartoon of the dissipative dynamics. The pumping dynamics can be visualized by dissipative 
pushing of excitations (green and red dots) between adjacent corners of the plaquette. 



